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invariant  for  a  large  class  of  optical  systems.  It  is  also  shown  that  fundamental 
limitations  for  the  concentration  of  light  follow  from  the  uncertainty  principle  and  the 
second  law  of  thermodynamics,  which  apply  quite  generally.  These  relations  show 
why  quasi-homogeneous  light  (such  as  light  from  thermal  sources)  can  not  be 
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1.  INTRODUCTION 


Non-imaging  concentrators  are  almost  always  studied  using  geometrical  optics,1 
because  no  wave  theory  techniques  have  been  developed  that  are  adequate  for  this 
purpose.  Recently,  some  progress  has  been  made  toward  developing  a  wave  theory  that 
might  be  useful.2,3,4  A  primary  difficulty  is  the  proper  definition  of  the  radiance 
function  within  a  wave  theory.  Walther5,6  has  suggested  two  different  definitions, 
which  have  received  some  study.7,8  In  this  report,  some  preliminary  observations  are 
made  on  the  way  such  a  theory  might  be  used  to  analyze  non-imaging  concentrators. 

The  concentration  of  light  of  arbitrary  coherence  is  considered  in  section  2.  An 
equation  is  derived  that  describes  the  propagation  of  the  radiance  function  for  such  a 
light  field  through  any  linear,  stationary  optical  system.  Some  limitations  for  the 
concentration  of  light,  which  are  usually  obtained  from  geometrical  optics,  are  derived 
using  this  wave  theory.  In  section  3,  the  special  case  of  light  from  a  quasi-homogeneous 
source  is  considered,  and  it  is  found  that  the  propagation  expression  for  the  radiance 
function  is  simplified.  For  a  class  of  linear,  stationary  optical  systems  (including  free 
space  propagation  in  a  short- wavelength  limit),  the  radiance  function  is  invariant  upon 
propagation  of  the  light  through  the  system. 


2.  RADIATION  IN  ANY  STATE  OF  COHERENCE 


We  will  begin  by  considering  a  light  field,  of  arbitrary  coherence,  propagating  away 
from  the  z  =  0  plane  to  the  right,  as  shown  in  Figure  2. 1 .  For  convenience,  the  energy 
transported  by  each  monochromatic  component  of  this  field  is  described  using 
the  complex  form  of  the  radiance  function  over  the  z  =  0  plane,  as  defined  by 
Walther  (W),6 


'W.  T.  Welford  and  R.  Winston,  High  Collection  Nonimaging  Optics,  Academic  Press,  New 
York  (1989). 

*W.  Welford  and  R.  Winston,  “Generalized  Radiance  and  Practical  Radiometry,”/.  Opt.  Soc. 
Am.  A:  Opt.  Image  Sci.  4,  545-547  (1987). 

JR.  Winston  andX.  Ning,  “Generalized  Radiance  of  Uniform  Lambertian  Sources,”/.  Opt.  Soc. 
Am.  A:  Opt.  Image  Sci.  4, 516-519  (1988). 

*R.  Winston,  “Analogy  Between  Fresnel  Diffraction  and  Generalized  Radiance,”/.  Opt.  Soc. 
Am.  A:  Opt.  Image  Sci.  6, 145-146  (1989). 

5A.  Walther,  “Radiometry  and  Coherence,"/.  Opt.  Soc.  Am.  58,  1256-1259  (1968). 

‘A.  Walther,  “Radiometry  and  Coherence,”/.  Opt.  Soc.  Am.  63, 1622-1623  (1973). 

7J.  T.  Foley  and  E.  Wolf,  “Radiometry  as  a  Short  W avelength  Limit  of  Statistical  Wave  Theory 
with  Globally  Incoherent  Sources,”  Opt.  Common.  55, 236-241  (1985). 

'J.  T.  Foley  and  E.  Wolf,  “Radiance  Functions  of  Partially  Coherent  Fields,”/.  Mod.  Opt.  38, 
2053-2068(1991). 
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oa  oo 

B§)(x'2,s)  =  ?~  j  jH^0)(x(,X2)exp[-ifa  (x;-X2)]d2x; ,  (2.1) 

—  90  — 90 

where  W^^xJ.x^)  is  the  cross-spectral  density  function9  for  a  typical  monochro¬ 
matic  component  of  the  field  at  frequency  w  between  two  points  on  the  z  =  0  plane 
defined  by  the  radius  vectors  xj  and  x’2  from  the  origin,  s  is  a  unit  vector  from  the 
origin,  0  is  the  angle  that  s  makes  with  the  +z  axis,  and  k  =  u/c.  The  radiance  function, 
as  usually  defined  in  radiometry ,  is  the  energy  radiated  from  a  unit  area  about  a  point 
x2  in  a  plane  into  a  unit  of  solid  angle  about  some  direction  s  from  the  origin.  In  a  wave 
theory  it  is  not  possible10  to  define  a  function  that  has  all  of  the  properties  of  the 
geometrical  optics  radiance  function.  There  is,  however,  more  than  one  definable 
function  that  approximates  the  geometrical  optics  radiance  function  in  some  useful 
respects.  The  definition  for  the  radiance  function  in  Equation  2.1  differs  from  that 
defined  earlier  by  Walther5  and  used  extensively  by  Marchand  and  Wolf11  in  then- 
development  of  the  statistical  theory  of  radiometry  and  radiative  transfer.  No  reason 
has  yet  been  discovered  to  prefer  one  definition  over  the  other,  except  its  usefulness  in 
a  particular  theory. 


*W.  H.  Carter,  “Coherence  Theory ,”  inThe  Optical  Society  of  America  Handbook  of  Optics,  M. 
Bass  (ed.),  McGraw-Hill,  New  York  (in  press). 

,0A.  T.  Friberg,  “On  the  Existence  of  a  Radiance  Function  for  Finite  Planar  Sources  of  Arbitrary 
States  of  Coherence,"  J.  Opt.  Soc.  Am.  69,  192-198  (1979). 

"E.  W.  Marchand  and  E.  Wolf,  “Radiometry  with  Sources  of  any  State  of  Coherence,”/.  Opt. 
Soc.  Am.  64,  1219-1226  (1974). 
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In  Equation  2.1  the  radiance  is  defined  as  a  function  of  the  cross-spectral  density 
function  over  the  z  =  0  plane,  which  can  be  defined  by9 

(xl .  *2 )  =  ^(0)(*i  ,o))[^(0)  (xj ,  w)]  ^  (2.2) 

where  (x', «)  is  a  random  variable  describing  the  complex  amplitude  of  a  typical 
monochromatic  component  of  the  field  with  radial  frequency  w  at  a  point  given  by  the 
radius  vector  x'  in  the  z  =  0  plane.  The  angle  brackets  denote  an  ensemble  average. 

If  we  assume  that  the  sources  of  the  field  are  to  the  left  of  the  z  =  0  plane,  as  shown 
in  Figure  2.1,  then  it  follows  that  the  amplitude  of  the  field  anywhere  to  the  right  of  this 
plane  can  be  represented  using  an  angular  spectrum  of  plane  waves  in  the  form 


4>{x,u)=  J  j  A^  (p,q)exp[ik(px+ qy + mz)]dp  dq ,  (2.3) 

m  =  -\jl-p2-q2 ,  if  p2  +q2  <\, 
m  =  ^p2  +  q2- 1,  ifp2+q2>l. 


where  the  complex  amplitude  of  the  plane  wave  propagating  in  the  direction  given  by 
the  unit  vector  p  =  (p.  q,  m)  is  the  Fourier  transform  of  the  field  amplitude  over 
the  z  =  0  plane;  that  is, 


oo  oo 

/li0)(P-9)  =  A  |  J ^(0)(x',w)  cxp[-ik{px' + qy'j\dx' dy' , 


(2.4) 


which  has  the  inverse 


*(0)(x',u)  =  j  J  A(°\p,q)  exp [ik{px' + qy')\dp dq .  (2.5) 


Substituting  from  Equation  2.5  into  Equation  2. 1 ,  we  have  a  very  simple  expression 
for  the  radiance  function  of  this  field, 

B^)(x',s)  =  ^<0)(x',«)]  A^0) (sx,sy )J exp( i*s •  x ') cos 0 ,  (2.6) 

where  the  components  of  the  unit  vector  are  s  =  (sx,  sy,  st).  From  Equation  2.6,  we  see 
that  Walther’s  radiance  function  for  a  light  field  of  arbitrary  coherence  is  simply  the 
ensemble  average  of  the  product  of  the  complex  conjugate  of  the  field  amplitude  over 
the  z  =  0  plane,  the  angular  spectrum  of  plane  waves  for  the  field,  a  Fourier  kernel 
between  the  conjugate  variables  s  and  x',  and  the  cosine  of  6,  the  angle  that  s  makes 
with  the  +z  axis.  It  is  clear  from  Equations  2.4, 2.5,  and  2.6  that  the  variables  s  and  x', 
on  which  Walther’s  radiance  function  depends,  are  Fourier  conjugate  variables. 
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Figure  2.2  Propagation  of  the  radiance  function  through  a  linear,  spatially  stationary 
optical  system. 


To  treat  non-imaging  concentrators,  we  must  consider  the  propagation  of  the 
radiance  function  itself  through  some  arbitrary  optical  system  from  some  input  plane 
z  =  0,  (shown  to  the  left  in  Figure  2.2)  to  an  output  plane  z  =  0  (shown  to  the  right  in 
Figure  2.2). 

The  radiance  function  on  the  input  plane  is  given  by  Equation  2.6,  using  primed 
coordinates: 


Z#*°>(x',s')  =  (f^(x\u)]*  AwfeSytycx  p(i*s'x')cos0'  (2.7) 

The  radiance  function  on  the  output  plane  is  given  by  the  same  equation  with 
unprimed  coordinates: 

B^)(x,s)  =  ^[^(l=0)(x,w)]  A^0)(sx,sy)Jexp(iisx)cosff.  (2.8) 


If  the  optical  system  is  linear  and  stationary,  then  the  field  amplitude  in  the  output 
plane  is 


•o  oo 

^,(*=°)(X)W)  _ J  j>'=° \x',o,)h(x-x’)dx'dy'  (2.9) 


as  a  function  of  the  field  in  the  input  plane  and  the  point-spread  function  h(x)  for  the 
optical  system.  By  taking  the  Fourier  transform  of  Equation  2.9,  using  the  convolution 
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theorem , 12  and  then  using  Equation  2.4  twice  to  represent  the  Fourier  transforms  of  the 
fields  over  the  z  =  0  and  z'  =  0  planes  as  angular  spectra,  we  get 

ACJ=0)(p,q)  =  A^\p,q)h{p,q) ,  (2.10) 

an  expression  for  the  angular  spectrum  of  plane  waves  for  the  field  in  the  output  space 
as  a  function  of  the  angular  spectrum  in  the  input  space  and  the  transfer  function  for  the 
optical  system.  The  tilde  indicates  a  two-dimensional  Fourier  transform,  so  the  transfer 
function  is  given  by 


*(/»•«)  =^r  j  \h{*)e*9[-ik(px+qy)\dxdy,  (2.n) 


the  Fourier  transform  of  the  point-spread  function. 

Substituting  from  Equations  2.9  and  2. 10  into  Equation  2.8  and  using  Equation  2.7, 
we  get 


Bw  0)(x,s)  =  -iexp(/*sx){[fl<f'  0) (x,s)exp(-iks- x)] *  h*  (x)J ti(sx , sy )  (2.12) 


as  an  expression  for  the  propagation  of  the  radiance  function  through  the  optical  system 
from  the  input  plane  to  the  output  plane.  Thus,  the  radiance  function  transforms  by 
convolution  on  its  spatial  coordinate  with  the  point-spread  function  for  this  system, 
and  also  by  multiplication  with  the  transfer  function. 

Linear  systems  are  especially  easy  to  treat  Depending  on  the  form  of  the  point- 
spread  function,  this  system  may  or  may  not  be  a  non-imaging  concentrator.  If  the 
point-spread  function  reasonably  approximates  a  Dirac  delta  function,  then  the  system 
images  the  field  over  the  z'  =  0  plane  on  the  z  =  0  plane.  If  the  system  represents  free 
space  propagation  from  the  z  =  0  plane  to  a  z  =  0  plane,  which  is  parallel  to  it  and  some 
fixed  distance  away,  then  Equation  2.9  is  simply  the  Rayleigh  diffraction  integral  of 
the  first  kind13  and  Equation  2. 12  holds  with  h(x),  given  by  the  well-known  Rayleigh 
kernel  from  this  integral.  A  large  class  of  non-imaging  concentrators  are  linear  and 
space  invariant,  and  therefore  can  be  treated  using  this  theory  to  transform  the  radiance 
function.  However,  some  non-imaging  concentrators  transform  the  radiance  function 
in  a  much  more  complicated  way.  In  the  most  general  case,  we  can  not  give  a  specific 
equation  equivalent  to  Equation  2.12  for  the  transformation  of  the  radiance  function; 
however,  we  can  identify  some  restrictions  on  the  transformed  radiance  function  for 
even  the  most  general  concentrator. 

Since  the  two  variables  s  and  x  in  Equation  2.8  (and  s'  and  x'  in  Eq.  2.7)  are  Fourier 
conjugate  variables,  the  radiance  function  can  not  depend  on  them  independently.  If  we 
define  the  standard  deviation  for  each  Cartesian  component  of  these  variables  by 
equations  such  as 


13J.  Goodman,  Introduction  to  Fourier  Optics,  McGraw-Hill,  New  York  (1968). 

,JW.  H.  Carter,  "Three  Different  Kinds  of  Fraunhofer  Approximations:  I.  Propagation  of  the 
Field  Amplitude,"  Radio  Sci.  23,  1085-1093  (1988). 
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J"  j"jr2^(0)(x,w)dr  dy 

(A*)2=-^ - 

|  ji(0)(x,u)dx  dy 

— oo  —  oo 

oo  oo 

|  Jp2Ai0)(P>4)dPd<? 

J  {Ai0)(P’<?)4p^ 


J  JxV^^fx.tojdx  dy 

—oo  — oo _ 

oo  oo 

J  J^0)(x, <*>)<&■  dy 

—  oo  — oo 

oo  oo 

I  jpAi0)(p><?H>d<? 

—  oo  — oo 
oo  oo 

J  Ja10)(p-pVp^ 


(2.13) 


then  the  radiance  function  is  constrained  by  the  well-known  uncertainty  law  for  Fourier 
conjugate  variables;  that  is, 


vx(z)  =  *xAp>— , 
vy(z)  =  AyA^>^-, 


where  the  phase  space  volume  function 


(2.14) 


V  ~vx(z)vy(z)  =  AxAyApAq>(\/2ir)2  (2.15) 

has  at  least  some  of  the  properties  of  the  etendue,  which  is  the  product  of  the  area  of 
a  light  beam,  normal  to  its  direction  of  propagation,  and  the  solid  angle  that  the  beam 
subtends  with  the  beam  axis,  as  defined  by  geometrical  optics  theories.1  Since  the 
maximum  value  for  Ap  and  Aq  is  unity  (p  is  a  unit  vector),  we  see  immediately  from 
Equation  2.14  that  the  smallest  possible  spot  a  concentrator  can  achieve  is  given  by 

Sc  =  Min(A*  Ay)  =  ^  j  .  (2.16) 

Equation  2.6  shows  that  each  monochromatic  component  of  the  light  can  be  concen¬ 
trated  into  a  spot  with  a  radius  that  is  actually  somewhat  smaller  than  a  single 
wavelength.  Such  a  small  spot,  however,  is  not  always  obtainable,  as  shown  in  the 
following  discussion. 

Consider  the  light  concentration  properties  of  a  general  optical  system.  All  we  know 
about  this  system  is  that  the  phase  space  volume  function  for  the  radiance  function,  as 
it  propagates  through  the  system,  can  not  be  decreased.  This  is  always  true  of  any 
optical  system .  This  law  is  usually  derived  for  the  etendue  in  geometrical  optics  theory, 
but  in  our  wave  theory  it  follows  from  the  second  law  of  thermodynamics.  The  phase 
space  volume  represents  an  uncertainty  associated  with  the  light  field,  so  the  uncer¬ 
tainty  associated  with  the  output  light  field  can  not  be  smaller  than  that  associated  with 
the  input  light  field; 
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y(out)  AxAyApAq  ‘  ' 

This  law  places  a  fundamental  restriction  on  what  any  concentrator  can  accomplish 
with  a  given  light  field. 

Let  the  input  light  field  form  a  spot  centered  on  the  origin  over  the  2'  =  0  plane,  with 
radii  along  the  x  and  y  axes  defined  by  the  standard  deviations  Ax'  and  Ay',  respec¬ 
tively.  Assume  that  this  spotoflight  radiates  into  the  2'  >  0  half-space  within  numerical 
apertures  projected  onto  the  xz  plane  and  yz  plane  and  defined  by  the  standard 
deviations  A p'  and  A q'.  Since  p  and  q'  are  the  x  and  y  components  of  the  unit  vector 
p',  normal  to  the  plane  waves  that  make  up  the  radiated  field,  it  follows  that  ifAp'=  Aq', 
then  the  cone  of  light  is  symmetric  about  the  +2'  axis  and  makes  the  angle  O'  with  the 
axis,  such  that 


sin0'  =  Ap'=  Aq’.  (2.18) 

In  general,  we  will  assume  that  A p'  ¥=  Aq',  so  that  the  light  cone  is  not  symmetrical 
about  the  2'  axis. 

Let  the  general  optical  system  act  as  a  light  concentrator,  faming  a  spot  of  light 
centered  on  the  origin  in  the  output  plane,  with  radii  defined  by  the  standard  deviations 
Ax  and  Ay,  and  radiating  into  the  2  >  0  half-space  within  numerical  apertures  defined 
by  the  standard  deviations  Ap  and  Aq.  We  define  the  concentration  of  the  light  in  the 
usual  way  as  the  area  of  the  light  spot  before  concentration,  divided  by  the  area  of  the 
light  spot  after  concentration: 1 


CA  A f'Ay'  (2.19) 

Ax  Ay 

Substituting  for  the  standard  deviations  defining  the  spot  size  in  the  input  and  output 
spaces  from  Equation  2.15  into  Equation  2.19,  we  get 


A  VCm)ApAq 
V'(ou,)  Ap’Aq’  ’ 


(2.20) 


Thus,  upon  substitution  from  Equation  2.17  into  Equation  2.20,  we  have 


c^_ApA£_ 
Ap'Aq’  ’ 


(2.21) 


a  very  general  constraint  on  the  ability  of  any  optical  system  to  concentrate  light,  which 
is  mandated  by  the  second  law  of  thermodynamics.  Equation  2.21  states  that  the 
concentration  of  the  light  field  can  not  exceed  the  increase  in  the  product  of  the 
numerical  apertures  for  the  projection  of  the  light  cone  onto  thexz  planes  and  yz  planes 
before  and  after  concentration. 

For  the  special  case  of  an  input  light  spot  that  is  symmetric  about  the  2'  axis,  so  that 
Ap’  =  Aq'  =  sin  O',  and  a  symmetrical  optical  system,  so  that  Ap  =  Aq  =  sin  0,  we  find 
from  Equation  2.21  that  the  concentration  is  limited  by 
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c< 


(2.22) 


where  6  and  6'  are  the  angles  that  the  light  cones  make  with  the  z  axis  in  the  output  and 
input  spaces,  respectively. 

The  maximum  concentration  is  obtained  over  the  z  =  0  output  plane  if  the  light  cone 
converging  onto  this  plane  from  the  aperture  of  the  concentrator  makes  a  half-angle  8 
that  approaches  90°.  Thus,  the  absolute  limit  on  the  concentration  that  can  be  achieved 
by  any  light  concentrator  is  given  by 


C< 


1 

sin2  6' 


(2.23) 


This  limitation  is  identical  to  a  very  well  known  limit  on  the  concentration  that  is 
usually  derived  by  use  of  geometrical  optics. 1  The  present  derivation  shows  that  this 
limit  follows  from  the  second  law  of  thermodynamics  and  is  therefore  fundamental  to 
any  optical  system.  Equation  2. 16  gives  the  smallest  possible  spot  of  light  that  can  ever 
be  formed  with  any  light  field.  Once  we  are  given  a  specific  light  field,  however,  with 
some  fixed  boundary  condition  over  the  z'  =  0  plane  for  the  field  amplitude.  Equation 
2.23  gives  the  maximum  amount  by  which  the  given  spot  can  be  reduced  without 
violating  the  second  law  of  thermodynamics  (and  thereby  decreasing  the  entropy  of  this 
particular  field). 


3.  RADIATION  FROM  QUASI-HOMOGENEOUS  SOURCES 


Next,  we  consider  the  concentration  of  radiation  from  quasi-homogeneous  sources 
rather  than  more  general  kinds  of  radiation.  Thermal  sources  and  many  other  naturally 
occurring  (nonlaser)  sources  are  quasi-homogeneous.  For  radiation  from  such  sources 
it  is  more  convenient  to  begin  with  Walther’s  earlier  definition  for  the  radiance 
function,  which  has  been  used  extensively  by  Marchand  and  Wolf  (MW)11  in  their 
theory  of  radiometry.  This  definition  has  been  shown  to  be  fully  equivalent  to  the 
definition  used  in  Equation  2. 1  for  quasi-homogeneous  light  fields.14  For  a  field  from 
a  quasi-homogeneous  source  over  the  z  =  0  plane,  we  define  the  complex  radiance 
function,  as  shown  in  Figure  1,  by 


o(°) 

°MW 


(x;,s)=^  j  j  w,<0)^x;-i-y,x;-yjexp(-ite  x:)d2x:. 


(3.1) 


Fora  quasi-homogeneous,  secondary  source  over  the  z  =  0  plane,  the  cross-spectral 
density  function  over  the  z  =  0  plane  is  given  by  the  product14 

l4W.  H.  Carter  and  E.  Wolf,  “Coherence  and  Radiometry  with  Quasihomogeneous  Planar 
Sources,"  J.  Opt.  Soc.  Am.  67, 785-796  (1977). 
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v^0>(*i.*2)=/(0)«)m(0)(x:). 


(3.2) 


where  /(0)(x+)  is  the  intensity  distribution  over  this  plane,  and  -Xj)  is  the 

complex  degree  of  spectral  coherence  between  two  points  x(  and  x2  on  the  plane  at 
frequency  u.  For  a  quasi-homogeneous  source,  we  require  that  the  intensity  /(0)(x'+) 
beaslow  function  relative  to  /i(0)(xl), meaning  that  ^t(0)(xl)  has  a  bounded  domain 
of  support  and  that  /(0>  (x  \ )  must  vary  so  slowly  as  to  be  essentially  constant  over  any 
domain  of  equivalent  area.  The  coordinates  on  the  opposite  sides  of  Equation  3.2  are 
related  by 


( x[  +  X2 ) 

(ri+tf) 

X+“  2  ’ 

<N 

II 

.  + 

II 

V 

1 

/  ,  X- 

Xl=X++Y> 

,  ,  xL 

.  y- 

X2  =  X+-y, 

y*-y+~- 

(3.3) 


Substituting  from  Equation  3.2  into  Equation  3.1,  we  get 

fiMw(x,s)  =  /(0)(x)ii(0)(s)cosfl  (3.4) 

for  the  radiance  of  the  radiation  from  a  quasi-homogeneous  source,  where  the  tilde 
again  represents  a  two-dimensional  Fourier  transform  of  the  form 


M<0)($)  =  ^2-  J  JV0)(x)exp(-i*i  x)rf2x. 


(3.5) 


Equation  3.4  is  equivalent  to  Equation  2.6  for  the  radiation  from  a  quasi-homogeneous 
source. 

Since  both  definitions  we  have  used  for  the  radiance  function  are  equivalent  for 
quasi-homogeneous  sources.  Equation  3 .4  is  also  true  for  Walther’s  second  definition 
of  radiance,  so  that 


B$(x,  s)  =  /(O)(x)ji(O>(s)cos0 ,  (3.6) 

as  shown  specifically  by  Carter  and  Wolf 14  and  as  used  extensively  since  then.7-8  Thus 
Equation  2.12,  which  was  derived  in  general  for  any  type  of  radiation,  can  be  used  to 
study  the  propagation  properties  of  the  radiance  function  for  quasi-homogeneous  light 
as  given  by  Equation  3.6.  Before  discussing  this  further,  however,  the  significance  of 
Equation  3.4  will  be  examined  within  the  classical  theory  of  radiometry. 

The  radiant  intensity  radiated  from  the  z  =  0  plane  in  the  direction  given  by  the  unit 
vector  s  into  the  far  field  was  found  by  Marchand  and  Wolf !  1  to  be  given  by 
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oo  ao 

/u(s)  =  cos0  J  Jfl(0)(x,s)d2 


(3.7) 


where  for  a  quasi-homogeneous  source  we  can  now  drop  the  subscript  that  denotes 
which  of  Walther’s  definitions  we  are  using.  Thus,  substituting  from  Equation  3.4  into 
Equation  3.7,  we  find  that 


Jjs)  =  X2  cos2  0/'<O)(O)  A(0)(s) ,  (3.8) 


where 


/<0>(0)  =  t4  j  j/<0>M^ 


(3.9) 


is  a  constant  proportional  to  the  total  intensity  integrated  over  the  z  =  0  plane. 
Substituting  from  Equation  3.8  into  Equation  3.4  for  £(0) ,  we  have 


/w(x)  JJs) 
X2/(O)(0)  cos0  ' 


(3.10) 


From  Equation  3.10,  we  see  that  this  radiance  function  factors  into  a  function  of  x, 
which  is  proportional  to  the  intensity  distribution  over  the  z  =  0  plane;  a  function  of  s, 
which  is  proportional  to  the  radiant  intensity  radiated  from  the  plane  in  the  direction 
s;  and  a  1/cos  0  obliquity  factor.  Although  Equations  3.10  and  2.6  each  factor  into  a 
function  of  s  and  a  function  of  x,  note  that  s  and  x  are  not  Fourier  conjugate  variables 
in  this  theory  for  radiation  from  a  quasi-homogeneous  source  as  they  are  for  a  more 
general  theory  in  which  the  quasi-homogeneous  approximation  is  not  used.  From 
Equation  3. 10,  we  see  that  both  of  Walther’  s  definitions  for  radiance  can,  for  the  special 
case  of  quasi-homogeneous  light,  be  interpreted  as  describing  the  flow  of  light  energy 
from  points  in  the  z  =  0  plane  into  specific  directions  toward  the  far  field.  This  is  in 
keeping  with  the  usual  classical  idea  of  radiance. 

To  develop  a  wave  theory  for  a  non-imaging  concentrator,  we  again  consider  the 
transformation  of  the  radiance  function  from  the  z'  =  0  plane  to  the  z  =  0  plane  by  an 
optical  system,  as  illustrated  in  Figure  2.2;  however,  this  time  we  will  assume  that  the 
light  in  the  input  z'  =  0  plane  is  quasi-homogeneous. 

If  the  optical  system  shown  in  Figure  2.2  is  a  general,  linear,  stationary  system,  then 
we  know  that  the  cross-spectral  density  function  over  the  z  =  0  plane  is  given  by15 


^i=0)K.*2)=tt}  ]d\] 


xA(x,-x;)h*(x2-xi) 


(3.11) 


l5W.  H.  Carter,  “Generalization  of  Hopkins’  Formula  to  a  Class  of  Sources  with  Arbitrary 
Coherence,”  J.  Opt.  Soc.  Am.  A:  Opt.  Image  Sci.  2,  164-166  (1985). 
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as  a  function  of  the  cross-spectral  density  function  over  the  input  z  =  0  plane. 
Substituting  from  Equation  3.2  into  Equation  3. 1 1,  we  find  that  this  system  transforms 
the  cross-spectral  density  function  for  the  radiation  from  a  quasi -homogeneous  source 
into 


^=0)(x„x2)=^- j  px;/<*'=0)(*;) 

*1  j  ^2x1m<z/=0)(x-),,[xi-x+- y) 

x/t*fx2-x;+yj.  (3.12) 


From  Equation  3. 12,  we  see  that  the  field  over  the  output  plane  is  not  generally  quasi- 
homogeneous. 

There  are  some  linear,  stationary  optical  systems  that  can  transform  a  quasi- 
homogeneous  field  over  the  z  -  0  input  plane  into  a  quasi-homogeneous  field  over  the 
z  =  0  output  plane.  Transforming  coordinates  in  the  kernel  of  the  propagation  integral 
in  Equation  3.12  using  Equation  3.3,  we  have 


Thus,  we  find  from  substitution  of  Equation  3.13  into  Equation  3.12  that  to  produce  a 
quasi-homogeneous  field  over  the  z  =  0  output  plane  the  spread  function  must  obey  the 
equation 


+ t)/I*(X+  ~t)  =  •  (3.14) 

Substituting  Equation  3. 14  into  Equation  3 . 1 2,  we  find  that  the  the  two  surface  integrals 
in  this  equation  can  be  separated  so  that  the  intensity  transforms  between  the  two  planes 
as 


/(r_0)(x+)  =  -^- j  j/(r'"0>(x;)k/(x+-x;)<fr;dy;  (315) 


and  the  complex  degree  of  spectral  coherence  transforms  as 

oa  oo 

/I=0)(x_)=^- J  Jtt(z/“0)(x:)kM(x_-x:)dx:rfy:,  (3.i6) 
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where  the  cross-spectral  density  function  for  the  quasi-homogeneous  field  in  the  z  =  0, 
output  plane  is  given  by 

^=0)(x, ,  x2 )  =  /(-0)(x+ )  M(l=0)(x_ ) .  (3.17) 

Hereafter,  a  linear,  stationary  optical  system  with  a  spread  function  that  satisfies 
Equations  3.14, 3.15,  and  3.16  will  be  referred  to  as  a  quasi-homogenrous  preserving 
(QHP)  system.  An  example  of  such  a  QHP  system  is  the  linear,  stationary  optical 
system  with  a  Gaussian  spread  function  given  by 


Mx)  =  exp[_lx|2/(2<*2)]-  (3.18) 

Substituting  from  Equation  3.18  into  Equation  3.14,  we  find  that 

*/(x) =  — ~2  exp(  ix l2/<^2 ) »  (3.19) 

x  a  '  ' 

and 

^(x)  =  ^rexp[-,x|2/(4ff2)]-  (3.20) 

Thus,  any  linear,  stationary  optical  system  with  a  Gaussian  spread  function  is  a  QHP 
system. 

We  can  derive  an  equation  equivalent  to  Equation  2.12  for  the  transfer  of  the 
radiance  function  of  the  radiation  from  a  quasi-homogeneous  source  through  a  QHP 
optical  system.  Substituting  from  Equations  3.15  and  3.16  into  Equation  3.4  for  the 
radiance  function  over  the  z=0  plane,  and  again  using  Equation  3.4  to  represent  the 
radiance  function  over  the  z  =  0  plane,  we  get 


bmw0)(x-  s)  =  [*M^0)  (x>  »)  *  *1  (x)]  Ms) .  (3.21 ) 

where  the  tilde  indicates  a  Fourier  transform: 

1  oo  oo 

Ms)  =  ^2  J  jMx)exp[-i*(s-x)]d2x.  (3  22) 


Equation  3.21  is  identical  to  Equation  2.12,  which  was  obtained  for  a  light  field  of 
arbitrary  coherence,  except  that  *M(s)  is  not  the  Fourier  transform  of  tyx),  and  the  two 
Fourier  kernels  are  missing  from  Equation  3.21. 

The  minimum  spot  size  a  QHP  concentrator  can  produce  with  radiation  from  any 
quasi-homogeneous  source  can  also  be  calculated.  Because  the  field  in  both  the  input 
and  output  planes  is  quasi-homogeneous,  we  have  in  each  of  these  planes  an  intensity 
function  that  is  slow  relative  to  the  complex  degree  of  spectral  coherence,  so  that 
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Ax+  » Ax_  ,  Ay+  » Ay_  ,  ^  23) 

Axl  »  Axl ,  Ay'  »  Ayl , 

where  Ax+  and  Ay+  are  the  standard  deviations  of  A* =  ^(x^)  in  respect  to  x+  and  y+; 
Ax_  and  Ay_  are  the  standard  deviations  of  /*(,=0)(x_)  in  respect  to  x_  and  y_; 
Ax'  and  Ay'  are  the  standard  deviations  of  /(I  =0)(x' )  in  respect  to  x'  and  y' ;  and 
Axl  and  Ayl  are  the  standard  deviations  of  //*=0>(xl)  in  respect  to  xl  and  yl.  From 
Equations  3.4  and  3.5,  we  see  that  s  is  a  Fourier  conjugate  variable  to  x  ,  and  s'  is  a 
Fourier  conjugate  variable  to  xl.  Thus,  from  the  uncertainty  principle  for  Fourier 
conjugate  functions,  we  have 


A*-Ajx>^-,  Ay_Asy>A 

2x  21  (3.24) 

where  Ajx  and  Asy  are  the  standard  deviations  of  the  function  (s)  relative  to  sx  and 

jy,  and  As'  and  A s’y  are  the  standard  deviations  of  the  function  /j<I=0)(s')  relative  to 
sx  and  s'y.  From  Equations  3.20  and  3.21,  we  have 


Ax*. »  Ax_  > 
Ay+  »  Ay_  >  • 
Ax;»Axl  > 
Ayl »  Ayl  >  - 


X 

2xArx  ’ 
X 

2xA sy  ’ 
X 

2  *Asx  ’ 
X 

2lCASy 


(3.25) 


Thus,  the  smallest  spot  that  can  ever  be  achieved  with  radiation  from  a  quasi 
homogeneous  source  is  limited  by  the  equation 


SQH  =  Min(Ax+Ay+)»^-j  (3.26) 

for  the  maximum  values  of  Asx  <=  Asy=  1.  Comparing  Equation  3.26  with  Equation 
2. 1 6,  we  see  that  radiation  from  a  quasi -homogeneous  source  (such  as  a  thermal  source) 
can  never  be  concentrated  as  much  as  light  from  some  non-quasi-homogeneous 
sources  (such  as  a  laser  source). 

AlthoughEquation3.26givesabasic  limitation  on  the  smallest  possible  spot  of  light 
that  can  be  formed  with  radiation  from  a  kernel  source,  the  second  law  of  thermo¬ 
dynamics  gives  another,  independent,  limitation  on  what  can  be  done  to  concentrate  a 
given  radiation  field  from  a  quasi-homogeneous  source.  The  concentration,  as  defined 
in  Equation  2. 19,  can  be  applied  to  radiation  from  a  quasi-homogeneous  source  to  ob- 
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tain  the  optimum  concentration  that  can  be  obtained  by  a  QHP  concentrator.  To  obtain 
this  limitation,  we  again  use  the  second  law  of  thermodynamics,  as  in  section  2: 


yW  a  Ax+Ay+As'xAs'y  ^  (3.27) 

y(out)  Ax+Ay+AsxAsy 

so  that  upon  substitution  from  Equation  3.27  into  Equation  2.19,  we  get 


A x+Ay+  As' As'  ’ 


(3.28) 


which  is  equivalent  to  Equation  2.21.  If  we  assume,  as  before,  that  the  optical  system 
and  the  illumination  are  symmetrical  about  the  z  axis  in  both  spaces,  so  that 
As'  =  As'  =  sin  0'  and  Asx  =  Asy  =  sin  0,  then  we  find  that 


C<, 


(  sin  0 


sint 


so  that  the  absolute  limit  on  the  concentration  is  given  again  by 


(3.29) 


CS- (3-30) 
sur  0 

for  the  case  where  sin  0  approaches  unity. 

As  with  the  more  general  case,  we  have  two  independent  limitations  on  the 
concentration.  Equation  3.26  gives  a  fundamental  (diffraction)  limitation  on  the 
concentration  of  radiation  from  all  quasi-homogeneous  sources,  which  is  that  it  can 
never  be  as  well  concentrated  as  some  other  radiation  fields.  Equations  3.28  and  3.30 
give  another  limitation,  which  is  that  radiation  from  a  given  quasi-homogeneous 
source  can  not  be  concentrated  more  than  the  ratio  of  numerical  apertures,  before  and 
after  concentration,  without  violating  the  second  law  of  thermodynamics. 

Since  both  of  Walther’s  definitions  for  radiance  are  equivalent  for  quasi-homoge¬ 
neous  light,  there  is  another  approach  for  studying  the  propagation  of  the  radiance 
function  for  quasi-homogeneous  light  through  the  linear,  stationary  optical  system 
shown  in  Figure  2.2.  Substituting  from  Equation  3.6  into  Equation  2.12,  we  find  that 
the  radiance  function  for  the  field  in  the  output  plane  of  any  linear,  stationary  optical 
system  is  given  by 


fl(,=0>(x,s)  =  -y  jt(O)(s)cos0exp(i*sx))i(.sx,.sy) 

X 

x  |[/(0>(x)exp(-r*s  ■  x)]  *  A*(x)J ,  (3.31) 


provided  that  the  light  over  the  input  plane  is  quasi-homogeneous.  Consider  the  factor 
in  braces  given  by 
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F(x>s)  =  ^T{[/(0)(x)exP(-'*sx)]*A*(x)J 

09  OO 

=  J  J Ii0){x')h*{x-x')exp{-iks  x')  d2  x' .  (3 .32) 


If  we  assume  that  the  spread  function  h(x)  has  a  domain  of  support  about  the  origin  that 
is  very  small  relative  to  any  domain  over  which  /°>(x)  varies  appreciably,  so  that 
/f°)(x)  is  a  slow  function  relative  to  h(x),  then  Equation  3.32  can  be  approximated  by 


F(x,s)  =  I(0)(x)h*  (j ,, ,  sy  )  exp(-/*s  •  x) ,  (3.33) 

where  the  tilde  denotes  the  Fourier  transform  defined  by  Equation  2.1 1.  Substituting 
from  Equation  3.33  back  into  Equation  3.30,  we  have 

fi(z=0)(x,s)  =  £<0)(s)/<0)(x)cos 0|/i(s)|2  ,  (3.34) 

which  describes  the  manner  in  which  the  radiance  function  propagates  through  any 
linear,  stationary  optical  system  for  which  the  spread  function  satisfies  the  approxima¬ 
tion  that  t°\x)  is  a  slow  function  relative  to  h(x),  leading  to  Equation  3.33. 

Foley  and  Wolf78  have  shown  that  if  the  optical  system  is  simply  propagation 
through  free  space  from  one  plane  to  another.  Equation  3.34  is  closely  approximated, 
in  a  short-wavelength  limit,  by 

B£=0)(x,s)  =  A<O)(s)/<O)(x)cos0 ,  (3.35) 

so  that  the  radiance  function  is  the  same  in  both  the  input  and  output  planes  and  is 
therefore  invariant  upon  propagation.  From  Equation  3.34,  it  is  clear  that  for  propaga¬ 
tion  of  the  radiance  function  through  an  optical  system  that  has  a  unimodular  spread 
function  h( s),  so  that 


\h(sf  =  \  (3.36) 

over  the  domain  where  1st  £  1 ,  Equation  3.35  still  holds.  Thus,  there  is  a  whole  class  of 
optical  systems  for  which  the  radiance  function  is  invariant. 

The  free  space  propagation  of  a  field  that  contains  no  evanescent  plane  waves  has 
a  transfer  function  given  by 

h(s)  =  expjit/l-73  -rJ  (r  - r')] .  (3.37) 

which  is  unimodular.  But  the  corresponding  spread  function  is  the  Rayleigh  kernel, 
which  does  not  have  bounded  support.  Thus,  the  approximation  that  was  used  to  obtain 
Equation  3.33  does  not  hold  for  this  case.  Nevertheless,  Foley  and  Wolf  have  shown 
that  Equation  3.35  does  hold  for  free  space  propagation,  provided  the  wavelength  is 
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sufficiently  small.  The  radiance  function  for  light  is  therefore  invariant  for  both  free 
space  propagation  and  propagation  through  a  large  class  of  optical  systems. 


4.  CONCLUSIONS 


This  report  has  considered  the  concentration  of  a  light  field  in  any  state  of  coherence 
by  any  optical  system.  The  special  case  of  the  concentration  of  radiation  from  quasi- 
homogeneous  sources,  which  is  applicable  to  radiation  from  thermal  sources  and  most 
other  naturally  occurring  sources  of  low  coherence,  has  also  been  considered. 

A  simple  equation  has  been  derived  for  propagating  the  radiance  function  through 
any  linear,  stationary  concentrator.  This  relation  should  be  useful  for  modeling  such 
concentrators  when  a  wave  theory  is  required.  A  somewhat  simpler  equation  for 
describing  the  propagation  of  the  radiation  from  a  quasi-homogeneous  source  through 
a  QHP  concentrator  has  also  been  found,  which  reproduces  a  quasi-homogeneous 
source  field  in  the  output  plane.  Finally,  an  approximate  expression  has  been  obtained 
for  the  propagation  of  the  radiance  function  for  quasi-homogeneous  light  through  any 
linear,  stationary  optical  system.  This  relation  shows  that  for  a  large  class  of  systems, 
including  free  space  propagation,  the  radiance  function  is  an  invariant. 

It  was  possible  to  derive  very  general  limitations  on  the  smallest  possible  spot  of 
light  that  can  be  achieved  by  any  light  concentrator.  For  monochromatic  components 
of  a  light  field,  the  smallest  possible  spot  is  given  by  Equation  2. 16.  This  is  simply  the 
well-known  diffraction  limit  for  focusing  any  monochromatic  wave  field.  For  radiation 
from  quasi-homogeneous  sources,  it  was  found  in  Equation  3.26  that  the  smallest 
possible  spot  must  be  much  bigger  than  the  diffraction  limiL 

Finally,  it  was  shown  that  the  well-known  limitation  on  the  concentration  that  can 
be  achieved  with  a  given  light  source  can  be  derived  for  any  concentrator  and  any  type 
of  light  field  from  a  wave  theory,  and  is  a  requirement  of  the  second  law  of  thermo¬ 
dynamics. 
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